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The optical properties of the quarter-filled single-band CDW systems have been reexamined in 
the model with the electron-phonon coupling related to the variations of electron site energies. It 
appears that the indirect, electron-mediated coupling between phase phonons and external electro- 
magnetic fields vanishes for symmetry reasons, at variance with the infrared selection rules used in 
the generally accepted microscopic theory. It is shown that the phase phonon modes and the electric 
fields couple directly, with the coupling constant proportional to the magnitude of the charge-density 
wave. The single-particle contributions to the optical conductivity tensor are determined for the 
ordered CDW state and the related weakly doped metallic state by means of the Bethe-Salpeter 
equations for elementary electron-hole excitations. It turns out that this gauge-invariant approach 
establishes a clear connection between the effective numbers of residual, thermally activated and 
bound charge carriers. Finally, the relation between these numbers and the activation energy of dc 
conductivity and the optical CDW gap scale is explained in the way consistent with the conductivity 
sum rules. 
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I. INTRODUCTION 

The electrodynamics of low-dimensional multiband 
systems in which strong correlations responsible for for- 
mation of locaHzed electronic states coexist with cor- 
relations responsible for charge-density-wave ordering 

\ (CDW) still attract great attention. The vanadium chain 
compound BaVSs is a typical example of such (mag- 
netic) metals with the commensurate CDW instability 
[l|, H, 0, 01 . In this system, basic electrodynamic proper- 
ties are related to one wide band, while magnetic anoma- 
lies are associated with more or less localized states in two 
narrow bands. The relations between the mid-infrared 
(MIR) energy scales, measured at temperatures below 
the critical CDW temperature, the activation energies of 
transport coefficients and the related effective numbers 
of conduction electrons are not entirely clear [5] . In or- 
der to explain these relations in simple physical terms, 
one needs, in the first place, a precise description of the 
CDW coherence effects in the quarter-filled wide band. 

\ This problem that is still not completely solved for the 
single-band case is in the focus of this work. 

; The previous theoretical work [1 H, 0, B i, [M [H!, 

' m, [H, [3 [S m, [13, [II on the electrodynamics of 
the single-band CDW systems followed two basic routes. 
First, the collective contributions have been studied to 
explain the low- and zero-temperature properties of the 
incommensurate CDW systems, with the emphasis on the 
non-linear conductivity regime [ll,[il. The main contro- 
versies characterizing early theoretical investigations are 
resolved using accurate symmetry analyses, with a partic- 
ular care devoted to the local field effects [Tj] . The prin- 
cipal physical problem in these analyses is that the cou- 
pling between infrared-active collective modes and exter- 
nal electromagnetic fields is shown as a simple function of 
the single-particle (interband) polarizability [sl. [tI. [l3l [l3| . 



which implies the indirect electron-mediated, rather than 
the direct, photon-phase phonon coupling. The high 
(MIR) frequency analyses [13, [3], on the other hand, 
were focussed on the precise description of the single- 
particle excitations, in particular on the rigorous treat- 
ment of the square-root singularity in the optical con- 
ductivity spectra at fku = 2A. It is shown that the con- 
troversies related to the high-frequency optical spectra 
usually reflect the incorrect treatment of the diamagnetic 
current in the transverse response approaches, and, con- 
sequently, can be easily avoided using the longitudinal 
response theory [l9| . 

In this article, we consider a quarter-filled single-band 
CDW system and a related tetramerized metallic sys- 
tem in which all non-retarded short-range interactions 
are neglected, and all coherence effects in the correlation 
functions associated with the retarded, phonon-mediated 
interations are treated exactly through the known analyt- 
ical form of the current vertices and the inverse effective- 
mass tensor. Using this exactly solvable model, we re- 
examine several textbook results [1, [1, [3, [1| for the dc 
and optical conductivity of commensurate CDW systems. 
The analysis is however limited by the use of several 
approximations. We consider a Holstein-like electron- 
phonon coupling in which the coupling constant is inde- 
pendent of the electron wave vectors. Being interested in 
the optical conductivity spectra obtained by the reflec- 
tivity measurements, which are characterized by a rel- 
atively poor resolution in the low-frequency part of the 
spectra, we approximate the total optical conductivity 
by the sum of the collective and single-particle contri- 
butions, with the local field effects included only in the 
collective term. In the numerical analysis of the single- 
particle term, we consider the simplest case in which two 
relaxation rates are assumed to be independent of fre- 
quency and temperature. Finally, the commensurability 
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(Umklapp) effects are present in the model through the 
shift of the phase phonon frequency to a finite value and 
through the commensurate photon-phonon coupling (see 
Sec. 5). 

The article is organized as follows. In Sec. 2, the band 
dispersions of the tetramerized CDW case are briefly dis- 
cussed. In Sees. 3 and 4.1, we consider the electron- hole 
excitations in the two-band version of the model in which 
the Umklapp effects are neglected. The emphasis is on 
the clear identification of the effectively free and effec- 
tively bound charge carriers in the underdamped regime 
(the interband damping energy is small in comparison 
with the CDW gap) , their relation to both the activation 
energy of dc conductivity and the MIR scale in the opti- 
cal conductivity. The brief discussion of the optical con- 
ductivity spectra in the overdamped regime, tentatively 
related to the hybridizations/interactions with localized 
electronic states, is given in Sec. 4.2. In Sec. 5, the in- 
frared selection rules for the collective modes in the com- 
mensurate CDW systems are discussed in some detail. It 
is shown that the (interband) electron-mediated photon- 
phase phonon coupling vanishes for symmetry reasons, 
and that the direct coupling is proportional to the mag- 
nitude of the charge-density wave. 



II. TETRAMERIZED CDW CASE 

A. Bloch energies 

The starting point of the present analysis is the QID 
tetramerized model in which the electron-phonon cou- 
pling is formulated in terms of the displacement vector 



u 



(1) 

where (/i^ = Q • — is the position of an atom 

and (/) is an arbitrary phase. Two Fourier components 
of interest, u^q and u,pq, are functions of two ordinary 
components of the wave vectors q± Q, as follows u^q = 

(l/%/2) [Uq+Q-|-Uq_Q] and U^q = (i/V2) [Uq+q-Uq.q] . 

WAi and M^i correspond to two similar lattice deforma- 
tions which phase is shifted by 7r/2 and which frequency 
is expected to be strongly renormalized for low enough 
temperatures and for the band filling close to the quarter 
filling (evidently, u^q — u^_q, with \ = A.ip being the 
phonon branch index). Q (0.57r/a, Qy) is the nesting 
vector of the Fermi surface of the quarter-filled (CDW) 
case, or the tetramerization vector in a general (metallic) 
case. The related variations of the electron site energies 
(see Fig. 1), which are assumed to be the main mecha- 
nism of the electron-phonon coupling here, are 
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-y= ^ e"''''' [wipq sin 0^ -I- UAq cos 0i] . (2) 
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FIG. 1; Variations of tlie electron site energies (— (5e, +5e, 
+(5e, — 5e) induced by the periodic (tetramerization) lattice 
distortion u^q ^q,o(wA) and = 7r/4 in Eq. ((T}. The solid 
and dashed arrows illustrate the Raman-active and infrared- 
active q « modes. The effective charges of the white and 
the black ions are -\-q and — q, respectively. Here (5e = \/2A 
and q — qc/^/2. 



g is the electron-phonon coupling constant. The di- 
rect inspection of Fig. 1 shows that the periodic lat- 
tice distortion associated with the replacement UAq — *■ 
Sq,o{uA) + UAq leads to the charge transfer of the mag- 
nitude qc (determined in Sec. 5) modulated by the wave 
vector Q. This lattice distortion makes UAq and M^q to 
be the long wavelength phonon modes which are, respec- 
tively, Raman and infrared active. The other phonon 
modes are taken into account later in Sec. 3, and are 
treated on the same footing with the static disorder. 

The CDW instability of the quarter-filled QID elec- 
tronic systems in the case in which the disorder is ab- 
sent and the electron-phonon coupling g is independent of 
the electron wave vectors is described by the well-known 
textbook expressions (see Eqs. (|A1|) - (IA10[) in Appendix). 
The commensurability effects enter in the theory through 
the dependent (Umklapp) contributions in the Bloch 
energies i?L(k). For example, in the ID limit {ti, — > 0), 
the four Bloch energies ElO^) are {L = C,C, V,V_ is the 
band index) 

Sc,c(k) = -v/2i2+2A2±2i(k), 
Ev,v{^) = V^tl + 2A2 T 2t{k), (3) 
t(k) = y/2AHl + tl cos2 2k^a + A"* cos2 2(/), 

with ta and ti, being the bond energies in the highly con- 
ducting direction and in the perpendicular direction, re- 
spectively, and A = g{uA) / is the magnitude of the 
order parameter. When A is small enough, the disper- 
sions of two lowest bands can be approximated by 



Ec,c{^) = e(k)T^e2(k)_£^(k)£,(k)-C/(k). (4) 

Here £c(k) — e(k) — — 2^^ cos fcaOa and ec(k) — 
|£(k + Q)| are two relevant bare dispersions, and 2£(k) = 
£c{k)+ecik), e{k) « e{k) + A'^ / {2s{k)) are useful abbre- 
viations. The Umklapp processes are represented in the 
dispersions ^ through 



U(k) 
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(5) 



with 6{k) — A exp{i(/)(k)}, (/)(k) = sgn{kj.)(j> and a — 1. 
The quarter-filled a = hmit of Eq. ^ with £2(k) w 



£^ (k) + corresponds to the usual Peierls model which 
electrodynamic properties will be studied in more detail 
in Sees. 3 and 4. 



B. Optical conductivity in pure CDW systems 

In the present context, the main effect of the Umklapp 
processes is to shift the renormalized frequency uj^pQ of 
the q = (7a;ej; = phase phonon mode from zero to a fi- 
nite value 5]. Together with other pinning mechanisms, 
this gives the collective contribution to the optical con- 
ductivity of the form (see Fig. 2(a)) 



(a) 



m* uj[uj + i7^q) 



(6) 



with TO* being the collective mode effective mass. In the 
pure case, where the scattering of electrons on both the 
disorder and the phonon modes is neglected, Lu^q and 7^q 
in Eq. ^ are the solutions of the usual functional integral 
approach [10|, or the solutions of the Dyson equation 
shown in Fig. 2(b) 0,11, Ql- The structure of to* and i 



is discussed in more detail in Sec. 5 and Appendix. 

The related single-particle contribution to the optical 
conductivity (Fig. 2(c)) is given, for example, by the 
gauge-invariant expression 



(7) 



obtained by means of the longitudinal response theory 
[20t (the left-hand side of Fig. 2(c), with q — qaHa)- 
Xi,i(qa,^) is the sum of the intra- and interband charge- 
charge correlation functions defined in Appendix, and 
aijia^uj) is the related polarizability. 



III. BETHE-SALPETER EQUATIONS 

We take the collective mode contribution ^ aside and 
continue with the examination of the (single-particle) 
electrodynamic properties of the a — Q two-band model 
of Eqs. (HI),©. The Bloch energies are 



Ec,c{^) = e-(k) T y^(l/4)e2^(k) + A2 



(8) 



with ecc(k) = £c(k) — ec(k). To make the analysis more 
general, we consider two regimes of the two-band model 
close to the quarter filling in which an important role is 
played by thermally activated electrons in both transport 
and optical properties. The first regime corresponds to 
the ordered CDW state considered above, while the sec- 
ond regime is associated with the metallic state, where 
the doubled Fermi wave vector is 2fcF ^ Qx = 0.57r/a. 
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FIG. 2: The collective (a) and single-particle (c) contribu- 
tions to the optical conductivity in pure regime of the present 
CDW model. The solid lines are the electron propagators and 
the zigzag lines are the phonon propagators. The open circles, 
the open squares and the full circles represent the charge ver- 
tices, the current vertices and the electron-phonon coupling 
constants, respectively. The dashed line is the external scalar 
field V'^''*(q, oj) and the wavy line is the related electric field 
(see Eq. (|13|l ). (b) The Dyson equation for the phonon Green 
functions. The phonon self-energy is given by the expression 
(|A10|) in Appendix. 



Relaxation processes on the (quasi)static 
disorder 



Having in mind the known structure of the current 
vertices related to the bands ([8]) [2l|, or, alternatively, 
calculating the long wavelength charge vertices ([^5]) . we 
can make the generalization of the optical conductiv- 
ity ([7]) to the case in which the (quasi)static disorder 
is present. This can be easily done by using the longitu- 
dinal equation-of- motion approach [20', '25| , because this 
approach treats the (intraband) charge continuity equa- 
tion in a natural way (the calculation is more compli- 
cated in the transverse response theory [1^). The result 
(Eqs. ^ and ^) is equal to Eq. © with the adiabatic 
term -q in the charge-charge correlation functions replaced 
by two phenomenological damping energies /iF™*'''^ and 

The longitudinal equation-of-motion approach is one 
example of the calculations which describe the elemen- 
tary electron-hole excitations in the effective single- 
particle multiband models in the gauge-invariant way. 
The Bethe-Salpeter approach, described below, is an- 
other example. In both cases, the elementary excitations 
are represented by the Matsubara electron-hole propaga- 
tor 

2)"' (k, k+ , , k' , r) - ^ I?f2^') (k, k+ , , k' , r) , 



where 

Pf4(k,k+,k;,k',r) = i^ 
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FIG. 3: (a) The Bethe-Salpeter equations, (b) The related 
self-consistent equations for the induced densities 5n^^ (k, oj). 
The diamond represents the electron-hole self-energy, (c) The 
generalization of Fig. 2(c). 



and k+ = k + q. Here describes the scattering of 
conduction electrons on the (quasi)static disorder. Using 
this approach, it is possible to treat the relaxation pro- 
cesses on impurities, on lattice imperfections, and even 
on soft phonons. The hybridization with other (uncorre- 
lated) bands can also be studied in this way. 



B. Relaxation processes on phonons 

In the presence of inelastic electron scattering pro- 
cesses on phonon modes, or other boson modes in gen- 
eral, one can use a similar longitudinal diagrammatic ap- 
proach based on the self-consistent Bethe-Salpeter equa- 
tions for 

$^^'(k;ia;„) = $^^'(k, k+, k^, k', ic^„, ic^„+) 

(i(jj„+ = \LOn -l-ifc'n). In this case H'^ in Eq. ^ represents 
the electron-phonon Hamiltonian (|A4[) with the phonon 
wave vector q' running over the entire Brillouin zone and 
the phonon branch index A running over all (clastic and 
inelastic) scattering channels. (k;i(jj„) is defined by 
the Fourier transform of T)^^ (k, k+, k^, k', r). 



V 



LL' , 



k, k+,k'^,k',ii/„) 



(10) 



(?i//3)^$^^'(k,k+,kV,k',i 



n+) 



For temperatures and dopings not too close to the crit- 
ical values characterizing the CDW metal-to-insulator 
phase transition, the interband electron-phonon coupling 
in H[ can be neglected (notice that beyond this ap- 
proximation one needs four band indices in the electron- 
hole propagator Q). Now, for G'^'^(k -I- q',k) w 

q',k)|'« 

take the 



(k- 



VW(2M^ x<75^^(k + q',k) and 
|GA(q')| : the Bethe-Salpeter equations [2 



form (see Fig. 3(a)) 
[ihvn + i^LL' (k,k+) + I]L(k, itj„) - SL/(k+,iw„+)] 
x$^^'(k;iw„) = i[gL(k,itj„) -gL'(k+,ia;„+)] 



'k.k' 



l_ ^ |GA(qO 



x$^^'(k + q';icj„ + ij/„0 



(11) 



in obvious notation (_EiL'(k, k+) = -^^(k) — £'i,'(k+) is 
the electron- hole-pair energy). In the intraband chan- 
nel, these equations transform into two coupled equa- 
tions, one is the charge continuity equation and the 
other is the transport equation. These two equations 
have to be solved self-consistently. For this purpose, 
we retain in Eqs. (jlip only the most singular scatter- 
ing processes by putting iujn+ — > -E'i'(k+)/fi. and ia;„ 
ELO^J/fi, in two single-particle self-energies, I]L(k, iw^) 
and I]i/(k+, iaj„+), and in the related vertex corrections 
on the right-hand side of the equations. In this way, the 
expression in the brackets on the left-hand side does not 
depend on iwn, and the equations can be rewritten in 
terms of the intraband {L — L') and interband {L ^ L') 
induced densities Sn'"^ (k, uj) defined by 



(5n^^(k,cj) 



5]il?"'(k,k+,kV,k',. 

k' 



Here 



eg^'^(k+,k)F^ 



^'^(k+,k)I/-*(q,c.). 



(12) 



£^L'L (k+,k; 



\E^{lo), (13) 



with Ea{u)) and V^'^^{q,u}) — {i/qa)Ea{uj) being, respec- 
tively, the external electric field and the related scalar 
potential (notice that for L = L' Eq. (fT3|) reduces to the 
expression eV'^^*{q,uj) — e{i/qa)Ea{uj) that is indepen- 
dent of k). Furthermore, ■^(k+,k) is the charge ver- 
tex, given by Eqs. (lA8|) and ([Ml), «/a'^(k) is the current 
vertex and €[ = qaSa- 

In Eq. (dll), X'^^'(k,k+,k^,k',w) is obtained by the 

analytical continuation of I?^^ (k, k-|_, k'^, k', iun), \Vn 
UJ + irj. The summation of Eqs. (jlip over iw„ makes the 
single-particle self-energy contributions and the related 
vertex corrections to appear on the same footing, result- 
ing in the electron-hole self-energy of the form 



/iir^^'(k,c.) 



(14) 



E E 

Aq' s=+l,- 



|GA(qO| 

N 



q')) 



/^b(q') + /Gsi?L'(k 



nu} + iri + Ell, (k, k + q') + sh^xc^' 
/•j(q') + /(si?L(k + q')) 



nw + ir] + Ell' (k + q', k) - sfi^Aq' 
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The corresponding generalization of U^^ (k, w) to the 
case with the interband electron-phonon couphng is 
straightforward; however, it is beyond the scope of the 
present work, /^'(q') and /(-^^(k)) = /^(k) in Eq. ^ 
are, respectively, the Bose-Einstein and Fcrmi-Dirac dis- 
tribution functions. Finally, notice that the electron 
scattering on the static disorder, described by the po- 
tential t/i(q'), is included here through the replacement 

c^Aq' ^ and (1 + 2/b(q'))|GA(q')|V^ l^i(q')l' m 
the A = boson branch [2Q] . 



C. Transport equations and optical conductivity 



The resulting self-consistent equation for the densities 
6n'"^ (k, Lo] represents the generalization of the transport 



plicitly for the (single-band) case in which A = and 
the relaxation rate T is presumably small but finite 
(Ini{77"(k, cj)} w r). These expressions will make the 
comparison of the present approach with the common 
CDW approaches in Sec. 5 straightforward. 

For A = and (\~ qxe-x, the electron-hole-pair energy 
is e{}L+) — £(k) w fi^o(q) — fiqxVxi^)^ '^ith tuliq) being 
independent of k for k w kp. Eq. psp can be written 
now in the textbook form [2^ 1 

(jj(5no(k, uj) — LUo{q)Sni(k, uj) — 0, 
(cj + ir)(5ni(k, iv) — cJo(q)<5no(k, oj) 

df{k) 



equation 



-Bl'l (k+,k) 

In this equation Sn'"^ (k, lu) is the contribution to 
proportional to J„ ^ (k) (the effective current 
vertices in Eqs. (HH) and are Ji'^(k) = J^^Ck) in the 
intraband channel and Ja~(k) — ?iajJa~(k)/£'ii(k, k) 
in the interband channel). Equation (|15p is solved con- 
sistently with the charge continuity equation and com- 
bined with the definition of the optical conductivity ten- 
sor. The result is the optical conductivity of the effective 
single-particle multiband models 



ucti vity 01 tne ( 

[Hill mi, HI 



-^Lt^.Vi?L'L(k+,k) 



/i(k)~/i,(k+) 



hio + nn^L' (k, io) + Ell' (k, k) - 



i^i^L^ (k,k+) 



with TILL = 1 in the intraband channel and n^L = 2 
in the interband channel. In the dynamical limit, the 
last term in the denominator can be safely neglected. 
At this point it should be recalled that the current ver- 
tices J^^ (k) can be shown in simple physical terms 
as: J^^(k) = ew^(k) in the intraband channel and 
Jq— (k) = (i/fi,)Pi'— (k)i?LL(k, k) in the interband chan- 
nel; ^^(k) is the electron group velocity and P^— (k) is 
the related interband dipole vertex. It should also be no- 
ticed that the coherence factor in the textbook expression 
for the interband contributions [5^, 'sj is slightly different 
from that in Eq. (fT6)) . as well as that the phenomenolog- 
ical extension of the textbook expression is found to lead 
to the non-physical in-gap spectra (see Fig. 1 in Ref. [13] )■ 



D. 



A = case 



Before turning to the detailed numerical analysis of 
the expression let us write Eqs. and (fTB|) ex- 



(-) 



de{k) 



ievl{k)Ex{u). (17) 



Introducing the boson field "(/"(qj which time derivative 
is proportional to the current induced by the external 
electric field. 



:^(q,c^), (18) 



one obtains the equation of motion 



The resulting optical conductivity is given by the inter- 
band term in Eq. (|16p . i.e., by the expressions 



;(a;-hir) -c^g(q) 



(k)r (16) 



(20) 



Here 2?° (q, w) is the bare propagator of the field V'lq, ^) , 
is the effective number of conduction electrons, 



eff.O 



m 
V 



2dm 

deik) 



(21) 



^^(k) = {l/h)ds{k)/dkx is the bare electron group ve- 
locity and ruaa = 2ft^/(/:aa^) is the mass parameter. 



IV. DISCUSSION 

A. Underdamped regime 

For comparison with the textbook expressions for the 
single-particle contributions to the optical conductivity, 
it is convenient to turn back to the simplest case discussed 
in Sec. 3.1. The electron-hole self-energy is approximated 
by its imaginary part and the imaginary part is assumed 
to be small when compared to the interband threshold 
energy; i.e., 77^^'(k,w) w ir^^'(k,w) and r^^'(k,w) w 
r^^' (r^^ = r^^"" and = r^''*", hereafter). In the 
underdamped regime (F'^*'''' <ti 2A), this approximation 
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FIG. 4: Main figure: The normalized effective number of 
conduction electrons {maa/m)Vorf^ as a function of the band 
filling for 2ta = 0.5 eV, 2ti = Q, T ^ Q and Ao = (dashed 
line) and 10 meV (solid line). The diamonds represent the 
dimensionless electron group velocity [2/(fa7r^)]ii^ (ftp). In- 
set: The inverse effective mass tensors 7^i^(k) and 755~(k) 
for Qx ~ 0.5n/a and Ao — 10 meV (solid and dotted lines, 
respectively). The bare-band vertex (A = 0, dashed line) is 
given for comparison. 



does not affect the spectrum in a critical manner, neither 
in the metallic regime (with 2kp not too close to Qx), 
nor in the CDW state at temperatures well below the 
critical temperature. In the metallic case, the effective 
numbers of conduction electrons in the resistivity, or in 
the Drude part of the optical conductivity, consist of two 
parts, the residual one, which remains finite at T — > 0, 
and the thermally activated part. In the tf, ^ Peierls 
case there are only thermally activated electrons. Finally, 
in the Peierls case with t^/A not too small there is a 
portion of the Fermi surface on which the CDW gap is 
not developed, leading again to a finite residual number 
of conduction electrons. 

To compare briefly the correct estimations of these 
numbers with the oversimplified "semiconducting" ex- 
pressions encountered in the literature [26], we apply 
Eq. pop now to the simplest metallic situation corre- 
sponding to Q = (0.57r/a, 0) and tt — > 0. In this re- 
spect, it should be noticed that the present analysis is 
not focussed on the single-particle properties of the usual 
strictly ID model, the physics of which is dominated by 
the singular nature of the forward scattering processes 
in the single-electron self-energy. Instead, we consider 
the electron-hole excitations in the typical QID model 
in which the QID character is incorporated in the afore- 
mentioned form of the electron-hole self-energy, where 
the forward scattering contributions are removed by the 
Ward identity cancellation of the q self-energy and 
vertex corrections in Eq. (|14p . The ^6 — > limit serves 
here only to replace most of the integrations over the 2D 
Brillouin zone by the ID integrations. The results ob- 
tained in this way can be compared term by term to the 
well-known analytical expressions. 



FIG. 5: Main figure: The temperature dependence of the 
normalized effective number of conduction electrons in the 
(imperfectly nested Qy = 0) CDW state for 2ta = 0.5 eV, 
2tb — (solid line), 1 meV (dashed line), 10 meV (dot-dashed 
line) and 15 meV (dot-dot-dashed line), 2fcF — Qx ~ 0.57r/a, 
To = 66.7 K and Ao — 10 meV. The diamonds represent the 
function exp{-A(r)/(fcBT)}. 



For example, the intraband conductivity along the 
highly conducting direction is proportional to the effec- 
tive number of conduction electrons n'iZ, as seen from 



^mtra/', ,\ 



m(uj + ir™tfa) 



(22) 



for a — X. njj^ includes the contributions of both bands, 
and it can be shown in two alternative ways, as a surface 
integral in the k space, or the related volume integral. 



dEL(k) 
^E^"W/^('^)- (23) 



V 



Lk(T 



In the first representation, it is given in terms of the elec- 
tron group velocity, while in the second representation it 
depends on the dimensionless inverse effective-mass ten- 
sor (k) = {m/h)dv^{\s)/dka. The dependence of 
on the electron doping 5 = 2ak-p /tt is shown in Fig. 4 for 
two choices of the distortion potential, Aq = and 10 
meV, and T ^ 0. The electron doping 5 = nVo {n- is the 
concentration of electrons and Vo is the primitive cell vol- 
ume of the Aq = lattice) and the dimensionless electron 
group velocity are also given for comparison. The latter 
represents the analytical solution of the tb ^ Q model. 

The metallic regime of the primary experimental in- 
terest corresponds to the case of a low electron/hole 
doping with respect to the quarter-filled band. In this 
case, the residual part in (shown in Fig. 4) is small 
when compared to its A = value (dashed line in the 
figure). Even for relatively low temperatures, the ther- 
mally activated part in becomes important, leading 
to pronounced thermal effects in The thermal ef- 

fects in the Drude part of the optical conductivity, or in 
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FIG. 6: The development of the normalized total optical 
conductivity with temperature in the underdamped regime: 
2ta = 0.5 eV, 2tb = 1 meV, To = 66.7 K, Ao = 10 meV and 



5 meV. 



Another important observation is that the optical con- 
ductivity model pB]) has a quite general form and is thus 
applicable to variety of low-dimensional systems. In this 
respect one should notice that the values of the effec- 
tive numbers are generally a direct consequence of 
the coherent scattering of conduction electrons on the 
crystal potential or on the periodically distorted lattice, 
irrespective of the band filling. These numbers represent 
the numbers of effectively free charge carriers, with the 
thermally activated part being usually very sensitive to 
temperature, providing the estimation of the activation 
energy A, or the charge-transfer gap in the general case. 
The corresponding numbers of bound charges, on the 
other hand, are the numbers of excited electrons across 
the CDW gap Q or across the charge-transfer gap. The 
latter numbers depend on temperature as well, as seen 
from the interband part of Eq. (|16p , which can be rewrit- 
ten in the form 



the resistivity, are further illustrated in Fig. 5 for the im- 
perfectly nested CDW case (corresponding to the elec- 
tron doping 5 = 0.5 in Fig. 4). The temperature de- 
pendence of the activation energy A is parametrized as 
A(T) = Ao(l - T/ToY, with (i = 1/2. nf^ is shown for 
four different values of tb and compared to the empirical 
expression nf^{T) = exp{-A(r)/(fcBr)}n°f^(To) (the 
residual part in n®^ corresponds to the zero-temperature 
interception value). The dashed and dot-dashed lines in 
the inset of figure represent the partial contributions of 
the lower and upper band for 4 = 0, n|^(C) and n°^(C) 
in Eq. (|23|) . respectively. The latter two have to be con- 
trasted to the numbers 



l^[l-/c(k)], ne = ^5]/c(k) 



(24) 



invoked from the theory of the ordinary semiconductors. 
The numbers (|24p are often taken as a good approxima- 
tion for rfJliC) and n':f^{Q. However, the disagreement 
between the expressions ()24|) and experimental observa- 
tion is typically of two orders of magnitude [2^ . and, con- 
sequently, the analyses based on the numbers ([2^ have 
to be taken with reservation. For the parameters used 
in Fig. 5, we obtain a factor of 25. This disagreement 
is easily understood on noting that the semiconducting 
approach, Eq. ([M)) . completely neglects the coherence ef- 
fects in the inverse effective-mass tensor. In other words, 
by replacing the vertices (k) by unity, when the num- 
bers n°Jl{C) and n°^(C) reduce to nh and ric indeed, this 
approach does not take carefully into account the fact 
that the thermally activated electrons are related only to 
the states in the thermal window around the Fermi level. 
Notice in this respect the singular behavior of the vertices 
7^3^ (k) at 2kx — 0.57r/a, shown in the inset of Fig. 4, 
which dominates the behavior of for 2fcF ~ 0.57r/a. 
On the contrary, the common approaches based on the 
surface-integral representation of usually give rise to 
the correct results. 



winter/, ,\ 



\L0 — I 



^"'(k)| 



/i(k)-/i,(k) 



(k,k)- 



(25) 



The temperature dependence of the total optical con- 
ductivity <P (w) = crj!f ''(w) + is calculated for 
several underdamped regimes. The damping energies 
^pmtra ^F"''""^ are assumed to be temperature in- 
dependent, for simplicity, and small in comparison with 
the scale 2A. The results for the imperfectly nested case 
are shown in Fig. 6 and are essentially the same as the 
results of the analogous nearly half-filled dimerized band 
[T9I [25I . Although it is physically clear that at tempera- 
tures T well below V'^^"' is much larger than Fjf*'"'^, the 
spectra shown in Fig. 6 are calculated for F™*'''* = rjT'^'^j 
for clarity. The shift of the MIR maximum in the spec- 
tra (fiii^MiR) with temperature follows roughly the en- 
ergy scale 2A, giving an alternative method of estimat- 
ing A [5]. The spectra are also characterized by the in- 
gap optical excitations, mainly in the energy range from 
^Wmin « 2A - fi-Fjj"*"' to 2A. The cut-off energy fej^in is 
associated with the zero conductivity interception of the 
line obtained by the linear extrapolation of the interband 
conductivity. The square-root singularity at fiw = 2A is 
absent here, and, according to Ref. [l^, the same shape 
{lo) is expected down to zero temperature. 



of a-'"*^'' 



B. Over damped regime 

A completely different situation appears in the over- 
damped regime where the damping energy fiF"''"'' be- 
comes comparable to 2A and where the real part of the 
electron-hole self-energy starts playing an important role. 
Typical physical situation corresponding to this regime 
is the CDW instability of the magnetic metals, where the 
wide band electrons exhibit the CDW instability, and are 
usually scattered by the strongly correlated electrons in 
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FIG. 7: Main figure: The optical conductivity as a function 
of damping energy ftFJ;"*'^'^ at T = 42 K. Dotted, dashed, long- 
dashed and dot-dashed lines correspond to ftr"'''^'^ = 20, 40, 60 
and 80 meV, respectively. The other parameters are the same 
as in Fig. 6. The solid line in the total optical conductivity 
at T = 67 K. Inset: The related interband contributions (the 
parameters are the same as in the main figure). 



one or more narrow bands. As the present discussion 
of the general expression (fTH]) neglects the real part of 
the self-energy, it turns out to be inappropriate for the 
detailed quantitative analysis of the overdamped regime, 
but it satisfactorily explains the interplay between two 
mutually competing energy scales, hT"^^°'^ and 2A. Fig. 7 
shows typical results. The pos ition of the MIR peak is 
at huiMiR ~ \/4A2 -I- (?irjj"*'='')2. In this case, fiwmin is a 
rather complicated function of the model parameters, but 
it is close to the value obtained in the narrow band limit 
of (7^^^°'^ (ui) , where all interband excitations are assumed 
to appear close to the threshold energy EqcO^, k) w 2 A. 
In this case, one obtains 



Re{<r('^)} 



{huj ~ 2A)2 -I- (?ir'nter)2 
^ ihLu + 2A)2 + (nrj;it"-)^ 



(26) 



As a result, in the overdamped regime /zu^jnin increases 
with increasing SF™*"'', for example, from the value 
hun^in « A at ftF|f"='' = 5A to fiw,„i„ « 2A at hP^^" = 
lOA. Knowing the experimental values of ^iwmir and 
^iWmin, we can again estimate A and fiF™*'^'' and compare 
A to the activation energy of transport coefficients. 

As mentioned in the introduction, a clear evidence of 
the overdamped regime is recently found in the vanadium 
chain compound BaVSa fsl • BaVSs is the multiband sys- 
tem with the (commensurate) CDW instability in which 
one weakly correlated wide band and two strongly cor- 
related narrow bands coexist at the Fermi level. Thanks 
to the extensive experimental activity in this system in 
recent years [H, 0, S 01 ; BaVSa seems to be a good candi- 
date for detailed investigation of the overdamped regime. 




+ 



(c) 




FIG. 8: (a) The collective contribution to the optical conduc- 
tivity in the common ID models with the forward-scattering 
short-range interactions, (b) The collective contribution in 
the present model, (c) The Hopfield-like processes neglected 
in the present approach. 



V. COMPARISON WITH THE COMMON CDW 
THEORIES 

In the low-energy examinations of the CDW systems 
it is common to use different strictly ID models with the 
non-retarded short-range interactions in which not only 
the forward-scattering but also the backward-scattering 
processes in the electron-hole self energy are neglected 
These approaches are usually based on the 
bosonization procedure. For example, for the simplest 
case with only the forward-scattering interaction (74 taken 
into account, the optical conductivity is given by the RPA 
diagrams in Fig. 8(a). The first diagram is nothing but 
the ideal intraband conductivity of Eq. (|20p . The interac- 
tion (74 leads to the renormalization of the electron group 
velocity in Eqs. ([^D]) and (PT|) in the way that 



while the electron dispersions in 9/(k)/9e(k), together 
with the related density of states, remain unaffected by 
(74. The boson field ip{q,uj) is foimd to satisfy the equa- 
tion of motion (|19p . while the optical conductivity is 
found to be given by Eqs. with u°(k) Up 

The present QID model differs from these strictly 
ID models in three important points. First, the 
electron-electron interactions in question are the re- 
tarded, phonon-mediated interactions d, H, second, 
the collective infrared-active mode is the phase phonon 
mode, which participates in the total conductivity spec- 
tral weight with usually less than 0.5 % (e.g., m* is of 
the order of 400m in K0.3M0O3 and, finally, the 

external electric fields are assumed to be well below the 
critical field required to give the non-linear conductivity 

In this section, we want to determine the structure of 
the photon-phase phonon coupling, and to compare the 
obtained results with both the textbook results obtained 
within the same QID model I, i, 0, i and the results 
of strictly ID models [1, We will derive here, for the 
first time, a clear selection rule for the infrared-active 
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(b) 



FIG. 9: (a) The direct photon-phonon coupling and (b) the 
indirect, electron-mediated photon-phonon coupling in pure 
CDW systems. 



(Fig. 9(a)) is given by the ground-state average of the 
coupling Hamiltonian 



^i?,(-q) J2 '5^i3;^(k,q)cLk.CLk. 

Aq 



LL'kcr 



xXc{ 



(30) 



(and Raman- active) phonon modes. On the qualitative 
level, the arguments arc the following. For the quarter- 
filled case shown in Fig. 1, the electric fields couple to 
the dipole moment of the phase phonon mode (the dipole 
moment is carried by the bound (condensed) electrons). 
For A comparable to the band width, this mode is the 
ordinary infrared-active phonon mode, characterized by 
the effective ion charge q — qc/V^, the reduced ion mass 
M and the concentration of ions I/Vq. On the other 
hand, for A small in comparison with the band width, 
the electrons in question are only weakly bound and the 
magnitude of the dipole is now proportional to e, which 
is q multiplied by the ratio between the amplitudes of 
the electron and ion oscillations. In the limit A — > 0, the 
dipole moment vanishes. 

To discuss this question quantitatively, we neglect 
again the Umklapp processes in the electron dispersions 
and focus the attention on the a = two-band model 
of Eq. (HI). In this case, all relevant vertex functions are 
simple functions of the auxiliary phase V'(k) defined by 



tan(p(k) 



2A 



(27) 



In the long wavelength limit, the electron- phonon vertices 
and the charge vertices of Eqs. (|A5|) and (|A8|) become 

■ y^(k+) + y>(k) ^ . 
— — sm w — sm (p(k), 



,i0(k) 



cosiy9(k), 



g^^(k+,k) = g§^{k+,k) 



■ ■ ¥^(k) - ^{k+) 

ism w 0{qa), 



gr(k+,k) « ie'^^'^) 



g^'^(k+,k) « q^{k+,k)^l, 



9^^(k+,k) 



h qaJ§-'~^{k) 
e £'cc(k+, k) 



(28) 
(29) 



According to Fig. 9, the photon-phonon coupling consists 
of two contributions, the direct coupling and the indi- 
rect, electron-mediated coupling. In the commensurate 
case in which the Hopfield-like processes [1^ of Fig. 8(c) 
are neglected, the linear direct coupling between external 
electromagnetic fields {Ex{—(l) is the electric component 
of the field polarized along the highly conducting direc- 
tion) and two q w phonon modes under consideration 



{sa = — 1, = 1 and the indices in the vertex ^(k, k) 
are A — ip, (f = A). Here UxXq is the amplitude of the 

electron oscillations and 5P^^^ (k, q) is the related dipole 
vertex. For the incommensurate case, the direct coupling 
is given in a similar way in terms of the current vertex 



(k,q) 



i|L,i-(k,k)|...,. 



(31) 



Since the phase phonon field can be written as Uxi^q = 
2i(ii^)<Pq [5|, we obtain again that the current vertex 
is proportional to the time derivative (/Sq, and the re- 
lated equation of motion is analogous to Eq. p^ . with 
e, m, n°^'° and ix'o(q) replaced by e, M, 1/Vb and cji^q. 
These parameters enter into the optical conductivity ex- 
pression through two adjustable parameters (m* and w^q 
in Eq. ([S]), for example). 

As a result, the direct photon-phonon coupling con- 
stant amounts to the replacement of c^^/ka'^ikcr in Eq. (|30p 



by 



rCLkcr) = 5L,c5L'fifc{k)- This result shows 



that the coupling function of the (Raman-active) mode 
UAc^^ 5'^*^(k, k), vanishes, while the coupling function of 
the phase phonon mode Utpq is finite and, as mentioned 
above, is proportional to the magnitude of the charge- 
density wave qc, which is given by 



'j^^^'^^'4^{^)fc{k). 



(32) 



Notice that the fact that the phonons u^q in Eq. (|30p cou- 
ple to external electromagnetic fields through the func- 
tion 5^*^(k, k), rather than f/^'-^(k, k), is a consequence of 
the fact that the variations V\{ri) in Eq. ^ are propor- 
tional to the derivative of VL\i with respect to r j , rather 
than to the displacement vector 

In the Hopficld-like processes the photon momentum 
is shared between an electron-hole pair and a infrared- 
active phonon in the way described by the dipole density 
operator 



^ JPi;;^(k + q'-q,k,q' 



"^L'k+q'— qcr 



CLko 



LL'q'kcr 



in Eq. ((30|) . with g\ -^(k, k) in the dipole vertex replaced 
by dA^O^ + q' ~ k). The sum over states in the phase 
phonon branch (q') is restricted to the region where the 
the vertex 5^ ^(k -I- q' — q, k) of Eq. I|A5|) is non-zero. 

The indirect, electron- mediated photon-phonon cou- 
pling shown in Fig. 9(b), on the other hand, is given 
by the interband correlation function xTT'^i'^^^) defined 
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in Appendix, or by Eq. (4.3) in Ref. According to 
Eqs. (UHl) and the product of vertices g^'^(k,k+) 

and (7^~(k+,k) in this correlation function is propor- 

tional to exp{— (/)(k)} J— (k), i.e., it is proportional to 
the bare electron group velocity f°(k). w°(k) is an odd 
function of and, consequently, the coupling function 
X™^°''(q, w) vanishes for both phonon branches. The 
background of this result is the fact that the electron- 
phonon coupling g is independent of the electron wave 
vectors. A more complicated form of the infrared selec- 
tion rules and the photon-phase phonon coupling func- 
tions can be expected in the models in which the electron- 
phonon coupling is related to the variations of the bond 
energies [6^ (namely, g depends now on the electron wave 
vectors, introducing a new channel in the electron-photon 
couphng). 

Finally, it is important to notice that the common text- 
book analyses 0,13, 11] explain the photon-phase phonon 
couphng in the way (see Eqs. (4.3) and (4.7) in Ref. Q) 
which is formally equivalent to the replacement of the 
bare electron group velocity in x'"^°'^(q, it'n) by v^{kF)- 
(That is, both the single-particle interband conductivity 
and the photon-phase phonon coupling are given in terms 
of the dimensionless function /(w) defined in Ref. 0].) 
Such correlation function is finite indeed and is of the 
same order as the expression (|32|) . However, when an- 
alyzing complicated CDW systems, this approximation 
has to be taken with reservation. 



VI. CONCLUSION 

In this article, we have derived the infrared selection 
rules for the single-particle and collective contributions to 
the optical conductivity in the tetramerized CDW system 
with the electron-phonon coupling related to the varia- 
tions of electron site energies. We have shown that the 
relations between the effective numbers of charge carriers 
and the activation energy of transport coefficients, or the 
MIR optical threshold energy, observed in the ordinary 
Peierls systems are naturally explained in terms of the 
CDW coherence effects. The same coherence effects are 
expected to be important in the correlated multiband 
systems with the CDW instability, as well. The quan- 
titative analysis of the latter systems is left for future 
investigations. 
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APPENDIX A: VERTEX FUNCTIONS AND 
CORRELATION FUNCTIONS IN THE 
TETRAMERIZED QID ELECTRONIC SYSTEM 

In order to fix the notation, here we rewrite the text- 
book expressions d, [E S @ for the exactly solvable 
tetramerized CDW case in which the electron-phonon 
coupling is independent of the electron wave vectors 
(Eq. ^ in the main text). 

The total Hamiltonian reads a,s H = Hq + H[ , with 
the bare Hamiltonian and the electron-phonon coupling 
Hamiltonian given, respectively, by 



Ho 



5]e(k)cLck. + ^EE 



2M 

kcr i=±lq«iO 

'ge' 

q^Oko- 



i^q+iqPq+'Q 



^ni,^ L V J V 



(Al) 



Here e(k) — —J2a "^ta cos kaUa is the bare electron dis- 
persion. Pq+iQ is the phonon field conjugate to Uq+jq 
and Wq_|_;q is the bare phonon frequency. 

The result of the phase transformation of H, which 
leads to the replacement UAq ^q,o(wyi) + ""Aq, is the 
effective tetramerized single-electron Hamiltonian 



- E<'(k)cLc/'k., 



(A2) 



ll'ka 



with the index / = — 1, 0, 1, 2, and with running over 
the first Brillouin zone of the tetramerized lattice. The 
matrix elements in the Hamiltonian (jA2p are i?o'(k) — 
e;(k) ~ £(k+ZQ), representing the dispersions of four ar- 
tificially tetramerized bands, and HQ'^^{\i) — Aexp{i0}, 
with A = g{uA) /^/N being the interband hybridization 
term. The diagonalization of (jA2p leads to 



(A3) 



iJS = ^i?L(k)c[k,C,.k. 



Lkcr 



The £'L(k) are four Bloch energies of the present prob- 
lem (Eqs. ^ in the main text), and they, together with 
the total energy of the electronic system, depend on the 
phase <j) in ilp'^^k). 

The observation (see the discussion of Fig. 1 in Sec. 2.1) 
that, for (ua) 0, the lattice vibrations associated with 
UAq and u^q are Raman and infrared active, respectively, 
means that there will be the direct photon-phonon cou- 
pling similar to the photon-phonon coupling in ordinary 
semiconductors. In principle, one can also expect the 
electron-mediated photon-phonon coupling, and, as dis- 
cussed in detail in Sec. 5, in the common microscopic ap- 
proaches, exactly this coupling is regarded as the dom- 
inant coupling 0, 0, It is given in terms of the 
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dimensionless electron-phonon coupling constants, 
gf^(k+,k) in 



LL'ko- 



(A4) 

and the electron-photon coupling constants (current ver- 
tices). 

The coupling constant g^^ ^ (k+ , k) in Eq. (|A4[) is given 
in the usual way 



9k 



'^(k+,k) = ^[e'^^c/k+qa,i')c^ka + i,i) 
I 

+e-'^>'U^+^{l + \,L')Ul{l,L)\ (A5) 



(k+ = k-|-q). Here 4)a = 4> and 4>^ = 0-|-7r/2 are two use- 
ful abbreviations. The U\^{l,L) are the transformation- 
matrix elements defined by 



with q = QaCa- 

Finally, after neglecting the relaxation processes in the 
electron-hole propagators associated with the disorder in 
the system, the renormalization of the phonon frequen- 
cies in (jAip is given in terms of the phonon self-energy 



(A9) 



where the correlation function XA,A(qjii'n) has the form 



-k-l-iqcr 



cL = Ef^ka,i)ci 



Lko-- 



(A6) 



The coupling of electrons to external long-wavelength 
scalar fields has the form which is similar to Eq. (jA4p . 



q LL'ka 

(A7) 

Here the ^^^(k+jk) are the dimensionless monopole- 
charge vertices, in the intraband channel equal to unity, 
and in the interband channel proportional to the inter- 
band current vertex ^(k), 



g^'^(k+,k) = J2Uk+^{l,L')U^{l,L) 



(A8) 



XA,A(q,ii^n) 



1 |.9r(k„k)|^ 

LL'ka 



ihVr, 



EL{k) 



(AlO) 



In the same approximation, the electron-mediated cou- 
pling of phonons to the external scalar fields is given 
by XA,i(qj i'^n)- This correlation function is obtained 
by replacing the vertices |g^^(k+,k)|^ in Eq. (jAlOp 
by ■^(k+,k) q^^ (k, k_|_). Similarly, the charge-charge 
correlation function Xi.i{<\i '^i^n) is given by Eq. (jA10|) 
with \q^^ (k, k4.)p being the vertices in question. For 
all three correlation functions we obtain Xi■j{^.^^)^ hj^ 
{1, A}, by the analytical continuation of Xi.j(q, iM„). 
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